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Abstract 

We consider a controlled Schrodinger equation with a dipolar and a 
polarizability term, used when the dipolar approximation is not valid. 
The control is the amplitude of the external electric field, it acts non 
linearly on the state. We extend in this infinite dimensional framework 
previous techniques used by Coron, Grigoriu, Letter and Turinici for 
stabilization in finite dimension. We consider a highly oscillating con- 
trol and prove the semi-global weak stabilization of the averaged 
system using a Lyapunov function introduced by Nersesyan. Then it 
is proved that the solutions of the Schrodinger equation and of the av- 
eraged equation stay close on every finite time horizon provided that 
the control is oscillating enough. Combining these two results, we get 
approximate controllability to the ground state for the polarizability 
system with explicit controls. Numerical simulations are presented to 
illustrate those theoretical results. 

Key words : Approximate controllability, Schrodinger equation, polariz- 
ability, oscillating controls, averaging, feedback stabilization, LaSalle invari- 
ance principle. 

1 Introduction 

1.1 Main result 

Following [R] we consider a quantum particle in a potential V{x) and an 
electric field of amplitude u(t). We assume that the dipolar approximation 
is not valid (see IDBA"*"! IDKABj ). Then, the particle is represented by its 
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wave function ■ilj{t,x) solution of the following Schrodinger equation 

idtij = {-A + V{x))^ + u{t)Qi{x)^ + u{tfQ2ix)ij, x£D, 
if^\dD = 0, 

with initial condition 

V'(0,x) = x£D, (1.2) 

where D C M™ is a bounded domain with smooth boundary. The functions 
V,Qi,Q2 G C°^{D,M) are given, Qi is the dipolar moment and Q2 the po- 
larizability moment. For the sake of simplicity, we denote by L^, Hq and 
respectively the usual Lebesgue and Sobolev spaces L^(Z),C), Hq{D,C) 
and H'^{D,C). The following well-posedness result holds (see |Ca) ) by ap- 
plication of the Banach fixed point theorem. 

Proposition 1.1. For any tp'^ £ HqCiH'^ and u G Lf^cii'^, +00), M), the sys- 
tem \1. 1\) ' [T^) has a unique weak solution ip £ C^{[0, +00), HqCiH'^). More- 
over, for all t > 0, Wipit, .)\\l2 = \\tl)^\\i2 and there exists C = C{Qi, Q2) > 
such that for any t > 0, 

\\ij{t,-)\\H2 < 11/11^26^/0 K-)l + l«MI'dr_ 

Let S := I?/' G L'^{D,C);\\ip\\i2 = l} and (•, •) be the usual scalar 
product on L^{D,C) 

if, 9) = I f{xW)dx, for /,5 G L\D,C). 
Jd 

We consider the operator (—A + V) with domain Hq n and denote by 
{\k)keN* the non decreasing sequence of its eigenvalues and by {4'k)keN* the 
associated eigenvectors in S. The family {(j)k)keN* is a Hilbert basis of L^. 
Our goal is to stabilize the ground state. As the global phase of the wave 
function is physically meaningless, our target set is 

C := {c4> ; c £ C and |c| = 1} , (1.3) 

where (p := 0i. 

Let J^o := {k > 2; {QkP, c^k) + o} and Jq := \k > 2;(Qi</),0fc) = o}. 
We assume that the following hypotheses hold. 

Hypotheses 1.1. 

i) VA; G Jo, {Q2<P, (pk) 7^ i.e. all coupling are realized either by Qi or Q2, 

ii) Card{jQ) < 00 i.e. only a finite number of coupling is missed by Qi, 
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iii) Ai — Afc 7^ Ap — Xq for k,p,q > 1 such that {1, k} ^ {p, q} and k ^ 1, 

Remark 1.1. The hypothesis is weaker than the one in |BN| (i.e. Jq = 0). 
As proved in |NH Section 3.4], we get that generically with respect to Qi 
and Q2 in C°°{D, R), the scalar products {Qi(j), 4>k) and {Q2<P, (pk) are ah non 
zero. The spectral assumption iii) does not hold in every physical situation. 
For example, it is not satisfied in ID if y = 0. However, it is proved in |NH 
Lemma 3.12] that if D is the rectangle [0, 1]", Hypothesis 1.1 iii) hold gener- 
ically with respect to V in the set g := {V £ C°°{D,R) ; V{xi,. . . = 

Vi{xi) + --- + Vn{Xn), with Vk e C°°([0,1],M)}. 

As in [CGLT| . we use a time-periodic oscillating control of the form 

u(t,V) :=a(V')+/3(V')sin(^^y (1.4) 

Following classical techniques (see e.g. |SVM] ) of dynamical systems in finite 
dimension let us introduce the averaged system 

idti'av = {- A + V{x))'lpav + a{'4)av)Ql{x)ll}av 

+ {a{i^av? + \l3{'<Pav?^ Q2{x)tPav, (1.5) 
'^av\Qj^ — 0, 

with initial condition 

V'a.(0,-) = V-"- (1-6) 

Let P be the orthogonal projection in onto the closure of Span {4>k \ k >2} 
and 7 be a positive constant (to be determined later). 

Our stabilization strategy relies on the following Lyapunov function (used in 
[HH]) defined on S* n Hi n i/^ by 

£(V) := 7ll(-A + V)Pnh + 1 - \{^^4>)?- (1.7) 
This leads to feedback laws given by 

a(V'a^,(i, •)) := -kh{iPav{t, ■)), (3{^av{t, •)) := gihi^avit, ■)), (1.8) 
with k > small enough and 

g G C^(M, M+) satisfying g{x) = if and only if x > 0, g' bounded, (1.9) 
and for j £ {1, 2}, for z £ H'^, 

Ijiz) = Im[7((-A + V)P QjZ, (-A + V)Pz) - {Qjz, <P) (</., z)\ . (1.10) 
We can now state the well-posedness of the averaged closed loop system 
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Proposition 1.2. Let R > 0. There exists = {V, Q2, R) > such that 
for any £ H"^ D D S with C{ip^) < R and k £ (0, ko), the closed-loop 
system lll.5\ }- n~S{ )- lll.8\} has a unique solution ij^av £ C''([0, +00), i?^ r]!!^). 
There exists M > such that 



Moreover, if AV'" G Hi n H"^ , then Aipav G C°([0, +00), n H'^). 

We define Xq := {^j^ £ S n n H"^ ; A?/;° e n H^} the set of admissible 
initial conditions. For an initial condition tp^ E Xq, we define the control 



where Va« is the solution of (fT3]l - (fL6]l - (fL8]l . 

The main result of this article is the following one. 

Theorem 1.1. Assume that Hypotheses {T71\ hold. LetC, the target set, be 
defined by U.'J\) . There exists ko = ^0(^)^2) > such that for any k € ko, 
for any s < 2 and for any t/j^ G Xq with < C{ip^) < 1, there exist an 
increasing time sequence (T„)„gN in M.*^_ tending to +00 and a decreasing se- 
quence (en)rieN in R*^_ such that if ip£ is the solution of U.l\ )- n~^) associated 
to the control defined by then for all n G N, if £ £ (0,e„), 



Remark 1.2. Theorem 11.11 gives the semi-global approximate controllability 
with explicit controls of system (|l.ip . Hypotheses 11.11 are needed to ensure 
that the invariant set coincides with the target set. The semi-global as- 
pect comes from the hypothesis < J~.{ip^) < 1 : by reducing 7 (in a way 
dependant of ip^), this condition can be fulfilled as soon as ip^ ^ C. 

In Theorem 11.11 there is a gap between the H^ regularity of the initial 
condition and the approximate controllabihty in H^ with s < 2. The extra- 
regularity is used in this article to prove an approximation property in H^ 
between the oscillating system and the averaged one (see Section [3|). Weak- 
ening this regularity assumption is an open problem for which an alternative 
strategy is required. The last lost of regularity comes from the application of 
a weak LaSalle principle instead of a strong one due to lack of compactness 
in infinite dimension. 

1.2 A review of previous results 

In this section, we recall previous results about quantum systems with bi- 
linear controls. The model (jl.ip of an infinite potential well was proposed by 
Rouchon in |Rj in the dipolar approximation (Q2 = 0). A classical negative 



\\i^av{t)\\H^ < M, Vt>0. 



(1.11) 




(1.12) 



dlstHs {Mt,-),C) <—, Vt G [Tn,Tn+i]. 
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result was obtained in |BMS) by Ball, Marsden and Slemrod for infinite di- 
mensional bilinear control systems. This result implies, for system (ll.ip with 
Q2 = 0, that the set of reachable states from any initial data in Pi H S 
with control in L^(0, T) has a dense complement in Pi Hq n S. However, 
exact controllability was proved in ID by Beauchard in [Bj for V = and 
Qi{x) = X in more regular spaces (iJ^). This resuh was then refined in |BL| 
by Beauchard and Laurent for more general Qi and a regularity H^. 

The question of stabilization is addressed in [BNj where Beauchard and 
Nersesyan extended previous results from Nersesyan |N1| . They proved, 
under appropriate assumptions on Qi, the semi- global weak stabilization 
of the wave function towards the ground state using explicit feedback control 
and Lyapunov techniques in infinite dimension. 

However sometimes, for example in the case of higher laser intensities, 
this model is not efficient (see e.g. [DB A"*" IDKABj ) and we need to add a 
polarizability term u{t)'^Q2{x)ilj in the model. This term, if not neglected, 
can also be helpful in mathematical proofs. Indeed the result of [BN] only 
holds if Qi couples the ground state to any other eigenstate and then the 
use of the polarizability enables us to weaken this assumption. Mathematical 
use of the expansion of the Hamiltonian beyond the dipolar approximation 
was used by Grigoriu, Lefter and Turinici in |GLT| [T] , A finite dimension 
approximation of this model was studied in |CGLT| by Coron, Grigoriu, 
Lefter and Turinici. The authors proposed discontinuous feedback laws and 
periodic highly oscillating feedback laws to stabilize the ground state. In this 
article, we extend in our infinite dimensional framework their idea of using 
(time-dependent) periodic feedback laws. We also refer to the book |Cor| by 
Coron for a comprehensive presentation of the feedback strategy and the use 
of time- varying feedback laws. 

How to adapt the Lyapunov or LaSalle strategy in an infinite dimensional 
framework is not clear because closed bounded sets are not compact so the 
trajectories may lack compactness in the considered topology. In this direc- 
tion we should cite some related works of Mirrahimi and Beauchard [BMI [M] 
where the idea was to prove approximate convergence results. In this article, 
we will use an adaptation of the LaSalle invariance principle for weak con- 
vergence which was used for example in |BN| by Beauchard and Nersesyan. 
There are other strategies to show a strong stabilization property. Coron 
and d'Andrea-Novel proved in |CdAN| the compactness of the trajectories 
by a direct method for a beam equation and thus the strong stabilization. 
In |CouH ICou2) Couchouron gave sufficient conditions to obtain the com- 
pactness in favorable cases where the control acts diagonally on the state. 
Another strategy to obtain strong results is to look for a strict Lyapunov 
function, which is an even trickier question, and was done for example in 
[GdANB] by Coron, d'Andrea-Novel and Bastin for a system of conservation 
laws. 

The question of approximate controllability has been addressed by var- 
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ious authors using various techniques. In |N2| . Nersesyan uses a Lyapunov 
strategy to obtain approximate contronabihty in large time in regular spaces. 
In |CMSB| . Chambrion, Mason, Si galotti and Boscain proved approximate 
controllability in for a wider class of systems using geometric control tools 
for the Galerkin approximations. The hypotheses needed were weakened in 
[BCCS| and the approximate controllability was extended to some spaces 



Explicit approximate controllability in large time has also been obtained 
by Ervedoza and Fuel in |EP) on a model of trapped ion, using different 
tools. 

1.3 Structure of this article 

As announced in Section [l.H we study the system (jl.ip by introducing 
a highly oscillating time-periodic control and the corresponding averaged 
system. Section [2] is devoted to the introduction of this averaged system and 
its weak stabilization using Lyapunov techniques and an adaptation of the 
LaSalle invariance principle in infinite dimension. 

In Section [3] we study the approximation property between the solution 
of the averaged system and the solution of (jl.ip with the same initial condi- 
tion. We prove that on every finite time interval these two solutions remain 
arbitrarily close provided that the control is oscillating enough. This is an 
extension of classical averaging results for finite dimension dynamical sys- 
tems. 

Finally gathering the stabilization result of Section [2] and the approxi- 
mation property of Section [3l we prove Theorem 11.11 in Section HI 

Section [5] is devoted to numerical simulations illustrating several aspects 
of Theorem 11.11 and of the averaging strategy. 

2 Stabilization of the averaged system 

2.1 Definition of the averaged system 

System (jl.ip with feedback law u defined by (|1.4p can be rewritten as 



where the operator A is defined by D{A) := n Hq, At/j := (iA — iV)il> 



in [Buq . 




(2.1) 



and 



F{s, z) := -i (a(z) + /3(z) sin(s)) Qiz - i {a{z) + (3{z) sm{s)f Q2Z. (2.2) 
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For any z, F{.,z) is T-periodic (with here T = 2ii). Following classical 
techniques of averaging, we introduce F^{z) := y Jq F{t,z)dt. We can 
define the averaged system associated to ()2.ip by 

f dttl^av = Alpav + F^{l{jav), 

Straightforward computations of F^ show that the system ()2.3p can be 
rewritten as (|1.5p . 

We show by Lyapunov techniques that we can choose a and /? such that 
the solution of the averaged system ()2.3p is weakly convergent in towards 
our target set C. 

2.2 Control Lyapunov function and damping feedback laws 

Our candidate for the Lyapunov function, C, is defined in (|1.7p . It is 
clear that C{ip) > whenever tp G S Cl Hq Ci and that J~.{tp) = if and 
only if -0 £ C 

The main advantage of this Lyapunov function is that it can be used to 
bound the norm. In fact, for any if; ^ S H Hq n H^, 

> 7ll(-A + V)Pi;\\l, > |||A(PV)lli2 -C> J||AV^||i, - C, 

where here, as in all this article, C is a positive constant possibly different 
each time it appears. This leads to the existence of C > satisfying 

M\ji2<cii + ciij)), y^JG snH^nH^ (2.4) 

Remark 2.1. Although the idea of using a feedback of the form (|1.4p is 
inspired by [CGLT] , the construction of the Lyapunov function and of the 
controls is here different because we are dealing with an infinite dimensional 
framework. We follow the strategy used in |NH iBNj . 

Choice of the feedbacks. ^We would like to choose the feedbacks a and 
/3 such that for all t > 0, —C{ipav{t)) ^ where ij^av is the solution of 

(USD, (112]). 

If I\il)av{t) G Hi n H"^ for alH > then 
^/:(V',„(t)) = 27Re[((-A + V)Pdt^av, (-A + 

- 2Re[{dt'4)av,4>){(t), i^av)] 

= 27Re ((-A + V)P{i^^av - iV^Jav - iaQitpav 

- i{o? + ^/32)Q2V'a.), (-A + y)PVa^ 



2Re 
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Then we perform integration by parts. As P commutes with (—A + V),V 
is real and thanks to the fohowing boundary conditions 

(-A + V)PiJav\dD = i^avldD = 4>\dD = 0, 



we have 



27Re[(-i(-A + VfPiJav, (-A + V)Pi;, 
-2Re^{{iA-iV)iJav,(l)){(l),i^c 
= 2jRe [(-iV(-A + y)PVa^, V(-A + 
+ 27Re ^{-iV{-A + V)PiPav, (- A + y)PVa.) 

+ 2AiRe {iipavA){(t>,^ai 
= 0. 



This leads to 

^£(V^a,(t)) = 2ah{^av{t)) + 2{a^ + i/?^^ h{^av{t)), (2.5) 
where Ij is defined in (|1.10p . 

In order to have a decreasing Lyapunov function we define the feedback laws 
a and /3 as in (jl.Sp . Thus (|2.5p becomes 

^^C{^l,a.{t)) = -2 (^Il{l - kh) - ^hg^h)^ ■ (2.6) 

If we assume that we can choose the constant k such that (1 — kl2) > 
for aU t > and if A^lJav{t) G Hq D H"^ then the feedbacks (|1.8|) in system 
(fT3]) lead to 

^^C{iPav{t)) < 0, yt> 0. (2.7) 

Well-posedness and boundedness proofs. Using the previous heuristic 
on the Lyapunov function, we can state and prove the well-posedness of the 
closed loop system (|1.5p - (jl.8p globally in time and derive a uniform bound 
on the norm of the solution. Namely, we prove Proposition 11.21 

Proof of Proposition I j.M By the explicit expression (jl.lOll of I2 , we get for 
any z G H'^, \l2iz)\ < /(|k||//2) where 

fix) := ||g2||L-+7(^+ll^l|L-+Al)(||Q2||c2X+||y||L-||Q2||L-+Al||Q2||L-). 
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Notice that / is increasing on M+. Let K := 2f yyj C{1 + Ryj where C is 

defined by ([231), := ^ and k £ (0, ko). 

The local existence and regularity is obtained by a classical fixed point 
argument : there exists T* > such that the closed loop system (ll.Sp 
with initial condition (|1.6p and feedback laws (|1.8p admits a unique solution 
defined on (0,T*) and satisfying either T* = +oo or T* < +oo and 

limsup\\ipav{t)\\H2 = +00. 

i-5-T* 

We have 

|/2(^a.(0))| < /(ll/lbO < / (^C'(1 + £(V/))) < 

thus, by continuity, \l2{tpav{t))\ ^ K for t small enough. 
Let 

Tmax ■ — sup {t G (0,T*);|/2(Va.(T))| <i^,Vr G (0,t)}. 
We want to prove tlicit TfYidx — 

T = +00. 

For all t G [0 

jTjYiax)-, we have (l — A;/2(^aii(^))) > 0, which implies (by (I2.6p ). 
C{^pavi')) is decreasing on [0,Tmax)- Estimate (|2.4p leads to 

llV'a.WII/f^ < + < yJc{l + C{ipo)), Vt G [o,r„„,). 

(2.8) 

Let us proceed by contradiction and assume that Tmax < T* . This implies 
\l2{,ilJaviTmax))\ = K. By definition of K, 

Ihii^avim < f (y^il + C^^ < ^ Vt G [0, Tmax). 

This is inconsistent with \l2{'4'av(Tmax))\ = K so T^ax = T* and the solution 
is bounded in when it is defined. As no blow-up is possible thanks to 
()2.8p we obtain that T^ax = T* = +oo and thus the solution is global in 
time and bounded. 

Finally, taking the time derivative of the equation we obtain the announced 
regularity. □ 

2.3 Convergence Analysis 

In all this section we assume that k G (0, A;o) where ko is defined in Propo- 
sition [T^] with R = 1. The closed-loop stabilization for the averaged system 
()L5p is given by the next statement. 

Theorem 2.1. Assume that Hypotheses li.il hold. If ip^ G Xq with < 
C{ip^) < 1, then the solution tpav of the closed-loop system il.5\} - [T78\) with 
initial condition U.6\) satisfies 

ipavit) ^ C inH"^. 
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We prove this theorem by adapting the LaSalle invariance principle to 
infinite dimension in the same spirit as in |BN| . This is done in two steps. 
First we prove that the invariant set, relatively to the closed- loop system 
(jl.5p - p.8p and the Lyapunov function C, is C. Here, Hvpothesis 11.11 is cru- 
cial. Then we prove that every adherent point for the weak topology 
of the solution of this closed-loop system is contained in C. This is due to 
the continuity of the propagator of the closed-loop system for the weak 
topology. 

2.3.1 Invariant set 

Proposition 2.1. Assume that Hvpotheses \l.l\ hold. Assume that ip^ belongs 
to S n Hq n and satisfies {tp^, cp) 7^ 0. // the function t 1— )• C{il)av{t)) is 
constant, then ip^ ^ C. 

Proof. Thanks to (g^]), the fact that (l - khi'ipavit))) > for alH > and 
([LQ]) we get 

hii^avi-)] = 0, l2{^Pav{-))9^ (^2 (V'a. (•))) = i.e. hii^avit)) > 0, Vt > 0. 

By p.Sp this implies that a{ipavi-)) = I^O'avi-)) = and then ipav is solution 
of the uncontrolled Schrodinger equation. So, 

oo 

Recall that (f> := (pi is the ground state. Following the idea of |N1) . we obtain 
after computations and gathering the terms with different exponential term 

where P{ip^ , j,k,Qi) and P{tp^,j,Qi) are constants. Then, by |NH Lemma 
3.10], 

{ij^,<P,){cP,tP^){Qi<P„4>)il+lX]) = 0, Vi G J^o- 

Using the assumption {(p,ip^) ^ and Hvpotheses 11.11 it comes that for all 
j g J-^Oi {'^^^4>j) = 0- This leads to 

V'a.(t) = e-^^^*(V°,</>><A+ J] e-'^^\,p\cp,)<p,, 
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where by Hypotheses ll.H Jq is a finite set. By simple computations we 
obtain, 

jeJo 

-K^°,</'>l'(Q2'/',0)) >0. (2.9) 

There exists Nq G N* and ipn)ne{o,...,No} such that 

n G {0,...,iVo}} = {±(Afe - A^) ; G Jq x (JoU{l})}, 

with ;Uo = and fij ^ if i 7^ k. Thus, (|2.9p imphes that for any n G 
{0, . . . , Nq}, there exists A„ = A„(^'^, (52) £ C such that 

^0 

Imf ^ Aje'''^M > 0, Vt > 0. (2.10) 

i=o 

Straightforward computations give 

^0 = E (7A|K</'„OP(Q2</'„ -KV'°,</')P(Q20,0). 

jeJo 

Thus, Im(Ao) = and our inequahty (|2.10p can be rewritten as 

No 

Im(^ J^^Aj-e^'^^*) > 0, Vt > 0, 
j=i 

with the //j being ah different and non-zero. Then using the same argument 
as in |CGLT| Proof of Theorem 3.1], we get that Aj = for j > 1 and then 
using ()2.9p in particular that the coefficient of g-^C-^j--**!)* vanishes. It implies 
{i/j^,(pj) = for all j G Jq- Consequently, ip^ = {ip^ ,(f))(f). As ip^^cj) ^ we 
obtain V'^ G C. □ 

2.3.2 Weak continuity of the propagator 

We denote by Ut{il)^) the propagator of the closed-loop system (jl.Sp - 
(jl.Sp . We detail here the continuity property of this propagator and of the 
feedback laws we need to apply the LaSalle invariance principle. 
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Proposition 2.2. Let Zn £ S n n H'^ be Qj seQuence such that — ^ Zqq 
in H^. For every T > 0, there exists N C (0, T) of zero Lebesgue measure 
verifying for all t G (0, T)\N, 

i) Ut{zn) Ut{zoo) in H"^ , 

ii) a{Ut{zn)) ^-^^a{Ut{zoo)) and j3{Ut{zn)) ^-^^ 13 {Ut{zoo)) ■ 

Proof. Proof of ii). We start by proving that if (z„)ngN £ FIq n satisfy 
Zn Zoo in H'^ tlien a{zn) — 5- a(-Zoo) and /3(zn) — >• (3{zoo)- Tlius ii) 

n— >oo n— >oo n— >oo 

will be a simple consequence of i). As proved in jBNl Proposition 2.2], using 
the fact that the regularity H^^'^ is sufficient to define the feedback, we get 

Ij{zn) — > Ij{zoo), for j = 1,2. 

So by the design of our feedback, 

aizn) — > a{zoo), l^izn) — > f3iZoo)- 

Proof of i). The exact same proof as in [BN, Proposition 2.2] based on 
extraction in less regular spaces, uniqueness property of the closed loop sys- 
tem and taking into account the polarizability term leads to the announced 
result. □ 

2.3.3 LaSalle invariance principle 

We now have all the needed tools to prove Theorem 12.11 

Proof of Theorem \2.1\ Consider G with < C{ip^) < 1. Thanks to 
the bound (|2.4p . lAt{ilJ^) is bounded in . Let (trj)neN be a sequence of 
times tending to +oo and -^oo G be such that ly(t„{tp^) ^oc in FI^- 

n— >oo 

We want to show that ipoo £ C. 

We prove that a(JAt{ipoo)) = and P(p(t{ipoo)) = 0. Indeed, the func- 
tion t ^ a{Ut{'>f)) belongs to L'^{Q,+oo) (because of ([22]) and ([L8])) so 
the sequence of functions [t G (0, -|-oo) i— )• a{lAt^j^t{i^^)) ^ tends to zero in 
-L^(0, -fcx)). Then by the Lebesgue reciprocal theorem there exists a subse- 
quence (trife)A;6N and A^i C (0, -|-cxd) of zero Lebesgue measure such that 

(/)) ^ 0, VtG (0,+oo)\iVi. 

fe— 5>00 

Let T G (0, -|-oo). Using Proposition 12.21 there exists N C (0, T) of zero 
Lebesgue measure such that 

a(Z^t+t„JV°)) a{Ut{^oo)), Vt G (0,r)\iV. 
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Hence, a{Ut{ilJoo)) = for all t G (0,r)\(iVi U N). The function t ^ 
a(Ut{ipoo)) being continuous we get a{Ut{i)oo)) = for all t G [0,T], and 
this for ah T > 0. Finally a(ZYt(V'oo)) = for alH > 0. 

The same argument holds for /3 as ^ : t i— )• l2{Ut{ilJ^)^ (I2 
belongs to -L^(0, +00). Then by the proof of Proposition 12.21 

~g{Ut+t^^{ij^)) ~9{Ut{^oo)), \fte{0,T)\N, 

and g[Ut{'ipoo)) = implies fi(Ut{ip 00)) = 0. 

These two results lead to the fact that C(lAt{4'oo)) is constant. 
By (pTl) . C{ipoo) < vC('0°) < 1 so {ipooA) / 0. All assumptions of Proposi- 
tion [271] are satisfied then ipoo G C. 

This concludes the proof of Theorem 12.11 and the convergence analysis of 
(P3|l . □ 



3 Approximation by averaging 

The method of averaging was mostly used for finite-dimensional dynam- 
ical systems (see e.g. |SVM| ). The concept of averaging in quantum control 
theory has already produced interesting results. For example, in [MSRJ the 
authors make important use of these averaging properties in finite dimension 
through what is called in quantum physics the rotating wave approximation. 
The main idea of using a highly oscillating control is that if it is oscillating 
enough the initial system behaves like the averaged system. We extend this 
concept in our infinite dimensional framework : we prove an approximation 
result on every finite time interval. More precisely we have the following 
result. 

Proposition 3.1. Let [s,L\ be a fixed interval and ip^ G Xq with < 
C{ip^) < 1. Let ipav be the solution of the closed loop system f fi.5]) . (775]) 
with initial condition ipavis, •) = . For any 6 > 0, there exists Eq > such 
that, if ip^ is the solution of associated to the same initial condition 

■0£(s, •) = and control u^{t) defined by with e G (0,eo) then 

\\A{t,-)-i^av{t,-)\\H^ <5, Vt G [S,L]. 

Remark 3.1. Notice that the controls a and (3 were defined using the averaged 
system in a feedback form but the control used for the system (jl.ip is 
explicit and is not defined as a feedback control. 

Remark 3.2. Due to the infinite dimensional framework, we are facing regu- 
larity issues and cannot adapt directly the strategy of |SVM| . 

Proof We define for (t, z, z) e R x H"^ x H'^ , 



F{t, z, z) := -i {a{z) + (3{z) sin(t)) Qiz - i (a(z) + /3{z) sm{t)f Q2Z. (3.1) 
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Notice that thanks to ([22]) for any {t, z) £Rx H^, 

F{t,z,z) = F{t,z). 



(3.2) 



With these notations the considered system (jl.ip with control ()1.12p and 
initial condition ip^^s,-) = can be rewritten as 



0, 



where ipav is the solution of the closed-loop system (jl.Sp with initial condition 
Denoting by Ta the semigroup generated by A, we have for any t > s, 



Mt)=TA{t-s)i;^+ I TA{t-T)F(-,^l^,{T),i;av{r)]dT, 



i^avit) = TA{t - S)^'" + / TA{t - T)F\iJav{T))dT. 

J s 

This implies for any t > s, 

\mt)-i^avmm < II fTA{t-T) [F(^,Va.(r)) -FO(^„,(r)) 



+ 



TA{t-T) F(-,V.(r),^„„(T) 



F(-,ijav{T) 



dr 
dr 



//2 

(3.3) 



We study separately the two terms of the right-hand side of (j3.3p . 

First step : We show the existence of C > such that for any t > s, for 
any e > 0, 



TA{t-T) F(^,Ve(r),Va.(r)) - F {J-,i;av{r) 



dr 



(3.4) 



<C/ W^eir) - ^av{T)\\fj2<^T. 
J s 

By (|2.2p . p.ip . it comes that for any r > s, for any e > 0, 

Ff^,Ve(r),V'a.(T)) -F(^,Va.(r)) 

= -i (a(V'a^,(r)) + /3(V'at,(r)) sin (J^^ Qi [^'^(r) - ^Jav{r)] 
i (a(V'a^,(T)) +/3(V'a^,(r)) sin (^))^(32 [V'£(t) - Va^r)] . 
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As ip av is bounded in ^ using (jl.Sp and (|1.10p we get the existence of 
Ml > such that for ah r > s, 



|a(^a.(r))| + \mav{r))\ < Mi. 



(3.5) 



As I sin (l^) I < 1, we get the existence of C > independent of e such that 
for any t > s, for any e > 0, 

|Fg,Ve(r),V'a.(r))-Fg,Va.(r)) \\^^ < C\\Mt) -^av{T)\\j,,. (3.6) 

Then the contraction property of imphes (|3.4p . 

Second step : We show that there exists C > satisfying for all t £ [s, L], 
for any e > 0, 



TA{t-T) ^F(-,Va.(r)) -FO(V',,(r)) 



dr 



< Ce. 



(3.7) 



We follow computations on the semigroup Ta done in |HL| . For {t,v) G 
M+ X C^{[s,L],H^), we define U and H by 

U{tM-))--= [\F{r,v{-))-F^{v{.)))dT, 
Jo 

H{t,v) : = dvU{t,v)v, 
where i) is the time derivative of v. 

Notice that the T-periodicity of F(-,v) and the definition of imply that 
U{-,v) is also T-periodic. 

Lemma 3.1. As ipav £ C^{[s, L], Hq n H'^), we have for any t G [s,L], for 
any e > 0, 



TA{t - t) [f ll,av{T)) - F\llja.{T)) 



dr 



eU {^-,^av{t)^ - eTA{t - s)U (^^,i^av{s)) 
+eA TA{t - t)U V'a.(T)) - e TA{t - t)H ^a.(r)) dr. 

Proof. The proof is done in |HLl Lemma 2.2] . 

We study separately each term of the previous right-hand side. 
• With K = [^J , we have ^ — kT G [0, T] and by periodicity 



□ 



(F(r,Va.(t))-i^°(^a.(t)))dr 

■t/e-KT 

F(r,Va.(t))-i^°(V'a.(i)))dr 
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As V' av is bounded in and a{ipav), /3(V'at)) s-r^ bounded there exists M2 > 
such that 

\\F{T,tPav{tmm < M2, \\F'^{iJav{t))\\m < ^2, Vr > 0,Vt > S. 

This leads to 



(■t/e-KT 

< I 2M2dr<2M2T, Vt>s,Ve>0. 







The same computations lead to 

TA{t - S)U (^-^,^av{s)) 



< 2M2T, Vt>s,Ve>0. 



Then, 



eU ( ^,iJav{t)] +eTA{t-s)U (J,Va^,(s) 



< Ce, Vt>s,Ve>0. 



(3. 



• By switching property, 

A TA{t - t)U V'a.(T)) dr = ^* rA(t - t)AU Va.(T)) dr, 
and for any t G [sj -^^]) for any e > 

AU (^^, Va.(t)) = Aj'^'' [F(r,V'a.(i)) - i^°(V'a.(t))] dr 

t/e~KT 

[AF{T,^av{t)) - AF^ijavit))] dr. 

By definition of F and F'^ we have 
AF{t, z) = -i{a{z) + I3{z) sm{t))A{Qiz) - i{a{z) + /3(z) sin(t))2^(Q2z), 
= -ia{z)A{Qiz) - i (^a{zf + \p{zf^ A{Q2z). 

By regularity hypothesis on Qi, Q2 and V there exists C > such that 

\\A{Q,z)\\h2 < C\\Az\\h^, \\A{Q2z)\\h2 < C\\Az\\h2. 

Thus thanks to Proposition 11.21 and the bound (j3.5p on a{ipav) and I3{ipav)i 
we get the existence of M3 > satisfying 



\\AF{T,^avm\H2 < M3, \\AF^iiJavm\H2 < M3, Vr > 0,Vt G 
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So, for any t G [s,L], for any e > 0, AU (^^jtpavit)) 
quently , there exists C > such that 



H2 



< 2M-^T. Conse- 



eA f TA{t-T)U (-,il;av{T))dT ^<Ce, Vt € [s, L], Ve > 0. (3.9) 
• For the last term we need to estimate H {^jipavit)) ■ We have 







(3.10) 



t/e-KT 

\dyF{T,'lpav)-dti'av " dF^ (l/jav) Alpav ) 

Using ([L8]) and (fTTOll . we have for any v,w £ C°([s, L], ifg^ n i^^ 
da{v).w = —kdli{v).w, d/3(u).u; = g' (^I2{v))dl2{v).w, 

where, 

d/j-(i;).u; = Im[7((-A + V)P{Qjw), (-A + F)Pu) 
+ 7((-A + V)PiQjv), i-A + V)Pw) 
- {QjW,4>){4>,v) - {QjV,(t)){(t),w 

Finally, we have 

dvF{t, v).w = —i{a{v) + /3(w) sin t)Qiw — i{da{v).w + d[i{v).w s\nt)Qiv 

—i{a{v) + (3{v) smt)'^Q2W — 2i{a{v) + f3{v) sin t){da{v).w + df3{v).w smt)Q2V, 

(3.11) 



and 



dF^{v).w = —ia{v)Qiw — ida{v).wQiv — i ( a{v)'^ + -/3(f )^ I Q2'w 



-i{2a{v)da{v) .w + /3{v)d/3{v).w smt)Q2V. 



(3.12) 



By Proposition O dt^pav G C'°([0, +oo),H^nH'^) so there exists M4 > 
such that 

\\dti;av(.t)\\H2 < Mi, VtG[s,L]. 

Hence the same computations as previously lead to the existence of C > 
satisfying 

\daii^av{t)).dti^av{t)\ + \dl3{i^av{t)).dti^avit)\ < C, Vi G [s,L], 

and thus by ()3.1ip . ()3.12p . for any t G [s,L], for any r > 0, 

\KF{T,i:av{t)).dti^av{t)\\m + WdF^i^avit)) Ai^avim m < C. 
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As a consequence, 



<CT, Vr G [s,L],Ve > 0, 



and then, 



e j\A{t - t)H (J-,^ljav{T)) \\^dT<{CLT)e. (3.13) 



We are now able to deal with the remaining term of the right-hand side of 
(j3.3p . Gathering inequalities p.Sp . (j3.9p and (j3.13p in Lemma [3.11 we obtain 
that there exists C > such that inequality ()3.7p holds. 

Third step : Putting together (j3.3p . (I3.4p and (j3.7p we obtain that there 
exists C > such that for any t G [s, L], for any e > 0, 

\\Mt)-i^avmH^<Ce + C f Ue{T) - i^av{T)\\H^dT. 



Hence Gronwall's lemma implies 

\mt)-^avmH^ < Cee^(*-^) < (Ce^(^-^))e, yt G [s,L], 

and Proposition 13 . II is proved with £q = — — r- . □ 

Remark 3.3. The proof we used is fundamentally based on the boundedness 
of Aipavit) on [s,L] and on Gronwall's lemma so it cannot be extended 
directly to an infinite time interval [s, +oo). 



4 Explicit approximate controllability 

The solution ipav of the averaged system (ll.5p . (|1.8p . can be driven in the 
weak topology to the target set C. The solution ip^ of the system (|l.ip 
associated to the same initial condition, with control w^, stays close to ipav 
on every finite time interval provided that the control is oscillating enough. 
Gathering these two results we prove Theorem ll.il 

Proof of Theorem \l.l[ We consider s < 2 fixed. 

By Theorem 12. 1| we can construct an increasing time sequence {Tn)nem 
tending to +oo such that for any n G N, 

distj,.(Va.(i),C) < Vt>T„. (4.1) 

Using Proposition 13 . 1 1 on the time interval [0,T„_|_i] we then construct a 
decreasing sequence {£n)n£N such that for any n G N, 

\\Mt)-i^avit)\\H' VtG [0,r„+i],VeG (0,£„). (4.2) 



5 NUMERICAL SIMULATIONS 
Then (gj]) , ([32]) imply that 

which is the statement of Theorem 11.11 

5 Numerical simulations 
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yte [r„,r„+i],Vee (o,e„), 

□ 



This section is dedicated to numerical simulations of system (jl.ip . First, we 
detail how we approximate the solutions of (ll.ip and (jl.Sp . Then, we check 
the validity of the implemented code. Finally, we illustrate different aspects 
of Theorem 11.11 and of the averaging property, Proposition 13.11 

5.1 Settings 

In all what follows, we set D = [0,1]. As the potential V will vary in 
this section, the eigenelements of —A + V are denoted fk,v ^ind Xk,v- Any 
function tp £ L'^{{0, 1), C) is approximated by its first M modes 

M 
k=l 

The unknown eigenvectors (pky ^re approximated in the following way 

N 

i=i 

The equality (—A + V)ipky = Xk.vy^ky leads to Ba^ = Xj^ ya^ with 

a'' = (ai , . . . , a%Y, B = diag(Ai,o, • • • , Xn,o) + {{V'fifl, ^j,o)) i<ij<N- 

Notice that Xj^Q = (fcTr)^ and ip^fl = -v/2sin(fe7r-) are explicit. The scalar 
products are approximated by the Matlab function quadl. The eigenele- 
ments and X^y are then approximated by the Matlab function eig. 

5.2 Approximation of ip^ and ipav 

Let 

i?o := diag(Ai,y, . . . , AM,y), Hn := {{Qnfiy,>fjy))^^ij^j^_^, nG{l,2}. 

It follows that the feedback laws (jl.lOp are approximated, for X G and 
iG{l,2}, by 

Ij{X) : = Im(7(Fo(0, {HjX)2, {HjX)MYy (^o(0, x^, • • • 
-{HjX)ixT), 
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leading to 

aiX) := -kh{X), /3(X) := - mm(/2(^), 0). 
Thus, if we define Xav, G M*^, systems (jl.ip and (jl.Sp are approximated 

by 

i^X, = {Ho + + ne(t)2F2)X„ (5.1) 

and 

i^^Xav = {Ho + a{Xav)Hi + {a\Xav) + ^P{Xa,f)H2)Xav, (5.2) 

where ite(t) = a{Xav{t)) + l3{Xav{t)) sin{t/e). Equations ()5.ip and (|5.2p are 
solved numerically (simultaneously) using Euler method with a time step dt 
and a Strang splitting method. 

5.3 Validation 

We now prove the validity of the implemented code. The eigenvectors ipky 
are approximated by = 50 modes. We take, as a test case, V{x) : = 
{x — 1/2)^, Qi{x) := x'^ and Q2{x) := x. The considered initial condition is 
■0*^ = + -^^2,v- The value of the oscillating parameter is e = 10~^. 

The parameter 7 is chosen such that C{^^) = 3/4. We compute the discrete 
Lyapunov function for the averaged system and the H'^ norm (with s = 1.8) 
to the ground state for both the oscillating and the averaged system. The 
time scale is [0,T] with T = 1000 and a time step dt = 10~^. For M = 5, 
we get the results presented in Figure [TJ 




Figure 1: Lyapunov function of the averaged system (left). H^ norm to 
the ground state (right) for the averaged system (continuous line) and the 
oscillating system (dashed line). 

As expected, we observe the convergence of the Lyapunov function to 0. The 
solutions of (jS.ip and (|5.2p are driven to the ground state (up to a global 
phase). To validate the simulations, we have also tested the code for M = 10 
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and M = 20. We obtained the same asymptotic behaviour and the same 
values for the Lyapunov function and the distance to the target. 
As the approximate controhabihty uses the fact that the controls are oscil- 
lating, the time step dt cannot be taken large with respect to the oscillating 
parameter e. For e = 10^^, we obtain the same results with dt = 10^^ 
and dt = 10~^. However, instabilities appear on the oscillating system for 
dt = 10~^. Thus, in all what follows the time step will be chosen smaller 
than e. We now present several simulations to illustrate various aspects of 
Theorem 11.11 

5.4 Influence of the initial condition 

For every other initial condition tested, the asymptotic behaviour is the 
same. We present here the results for the same parameters as in Figure [1] 
but with the initial condition ip^ = -^(fiy + -^^2,v + '~^^3,V- In this case, 
the stabilization of the averaged system is slower and we computed it for 
T = 5000. 




Figure 2: Lyapunov function of the averaged system (left). norm to 
the ground state (right) for the averaged system (continuous line) and the 
oscillating system (dashed line). 

We observe the same asymptotic behaviour as in Figure [TJ 

5.5 Averaging strategy 

We present numerically the influence of the oscillating parameter e. First, 
we consider the same potential, dipolar and polarizability moments as in 
Figure [TJ We compute the discrete norm (for s = 1.8) to the ground 
state (up to a global phase) and the discrete norm of Xav — Xf,. Figure 
[3] is obtained with e = 10~^ while Figure [His obtained with e = 10"'^. Both 
are computed with a time step dt = e and final time T = 500. 
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Figure 3: norm to the ground state (left) for the averaged system (contin- 
uous hne) and the oscillating system (dashed line) . gap from the average 
(right). 





Figure 4: norm to the ground state (left) for the averaged system (contin- 
uous line) and the oscillating system (dashed line) . gap from the average 
(right). 

For a fixed parameter e, we observe that the distance between the 
solution of (jl.ip and the solution of (jl.Sp with the same initial condition does 
not increase as the time goes to infinity but rather tends to a limit value. 
This limit value is of the same order of magnitude as e. We observe that 



\XayiT) - X,oMT)\\ 



\Xav{T) - Xio-4(r)||j^2 



30. 



This validates numerically the results of Proposition 13.11 and indicates that 
this averaging property should be valid on an infinite time horizon. 

The same behaviour has been obtained with other parameters. We 
present here the simulations with Qi{x) := cos{x) and Q2{x) := cos{2x), 
inspired by the physical situation of alignment dynamic of a HCN molecule 
as in |DKAB| . Figure [5] is obtained with e = 10~^ while Figure [6] is obtained 
with e = 10"^. Both are computed with a time step dt = e and final time 
T = 1000 (as the stabilization process seems slower in this case). 
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Figure 5: norm to the ground state (left) for the averaged system (contin- 
uous hue) and the oscillating system (dashed line). gap from the average 
(right). 





Figure 6: norm to the ground state (left) for the averaged system (contin- 
uous line) and the oscillating system (dashed line) . gap from the average 
(right). 



Remark 5.1. Although the time scales at stake in these simulations can seem 
very large, one has to remember that the Schrodinger equation is considered 
in atomic unity. 



6 Conclusion, open problems and perspectives 

In this article we have defined explicit oscillating controls that drive the 
solution of our system arbitrarily close to the ground state provided that the 
control is oscillating enough and the time is large enough. To achieve this 
we have used and developed tools from the theory of finite dimension dy- 
namical systems and applied them to the considered Schrodinger equation. 
We managed by adding a mathematically and physically meaningful term 
to weaken the previous assumptions on the coupling realized by this model. 
The assumptions that were made are proved to be generic with respect to the 
functions determining the system (potential, dipolar and polarizability mo- 
ments). The results presented should be generalizable to a compact manifold 



REFERENCES 



24 



with the Laplace-Beltrami operator. We performed numerical simulations to 
illustrate the approximate controllability. This gives numerical bounds on 
the time scale and on the values of the oscillating parameter needed to drive 
any initial condition arbitrarily close to the ground state. 

A challenging question would be to prove an approximation property of 
the averaged system on an infinite time interval [s,+oo). This would lead 
to approximate stabilization to the ground state. Based on the numerical 
simulations, this result seems to hold. Unfortunately the tools developed here 
are really based on the finite time interval and cannot be extended directly. 
In |BLj . Beauchard and Laurent proved the local exact controllability in 
around the ground state for the system (jl.ip in the dipolar approximation 
(i.e. Q2 = 0) under some coupling assumptions in one dimension. If one 
manages to extend their result to the system (jl.ip with suitable assumptions 
on Q2, this may lead to a global exact controllability result around the ground 
state, at least for the one dimensional case. The main difficulty would be 
to obtain the approximate convergence in the same functional setting as 
their local exact controllability result and with coherent assumptions on the 
polarizability moment. 

Acknowledgements : The author would like to thank K. Beauchard for 
having interested him in this problem and for fruitful discussions. 
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